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(57) ABSTRACT 

A global positioning system (GPS) receiver and method of 
acquiring and tracking GPS signals comprises an antenna 
adapted to receive GPS signals; an analog radio frequency 
device operatively connected to the antenna and adapted to 
convert the GPS signals from an analog format to a digital 
format; a plurality of GPS signal tracking correlators opera- 
tively connected to the analog RF device; a GPS signal acqui- 
sition component operatively connected to the analog RF 
device and the plurality of GPS signal tracking correlators, 
wherein the GPS signal acquisition component is adapted to 
calculate a maximum vector on a databit correlation grid; and 
a microprocessor operatively connected to the plurality of 
GPS signal tracking correlators and the GPS signal acquisi- 
tion component, wherein the microprocessor is adapted to 
compare the maximum vector with a predetermined correla- 
tion threshold to allow the GPS signal to be fully acquired and 
tracked. 
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RADIATION -HARDENED FAST ingly, there remains a need for a new, fully space-qualified 

ACQUISITION/WEAK SIGNAL TRACKING weak signal GPS receiver targeted for high altitude applica- 

SYSTEM AND METHOD tions. 


CROSS-REFERENCE TO RELATED 
APPLICATION 

This application claims the benefit of U.S. Provisional 
Patent Application No. 60/612,396 filed on Sep. 20, 2004, 
which is herein incorporated by reference in its entirety. 

GOVERNMENT INTEREST 

The embodiments of the invention described herein were 
made by employees of the United States Government, and 
may be manufactured and used by or for the United States 
Government for governmental purposes without payment of 
any royalties thereon or therefor. 

BACKGROUND 

1. Field of the Invention 

The embodiments of the invention generally relate to sig- 
nal processing and, more particularly, to techniques for signal 
acquisition and tracking. 

2. Background 

Global positioning systems (GPS) have found wide appli- 
cation for precision spacecraft navigation and formation fly- 
ing applications in Low Earth Orbits (LEO), but recent 
advances in GPS receiver designs and signal processing capa- 
bilities now make it feasible to consider using GPS to provide 
autonomous, onboard navigation capabilities for geostation- 
ary (GEO) or other high altitude space missions. 

Conventional (GPS) receivers are generally not readily 
capable of operating effectively at altitudes above LEO. This 
is primarily due, in part, to the fact that the GPS signals 
available there are approximately 10-100 times weaker much 
more sparsely present than on Earth or at LEO. However, GPS 
receivers typically provide a very accurate and economical 
means of navigation and have become very popular for LEO 
missions and could be just as popular for High Earth Orbit 
(HEO-above LEO) missions if the technology were available. 

Generally, GPS signal processing can be divided into two 
major tasks: signal acquisition (or detection) and then signal 
tracking (and demodulation). Signal acquisition is the pro- 
cess of finding a GPS signal within a two-dimensional 
unknown parameter space. Signal tracking is the process of 
continuously updating the estimates of these two and other 
signal parameters. Signal acquisition is generally a more dif- 
ficult process than signal tracking, and it is the inadequacy of 
conventional GPS acquisition techniques that generally pre- 
vents the use of weak signals in GPS receivers. 

The generally inadequate (for weak signals) conventional 
approach is to employ a serial search of the two-dimensional 
parameter space during acquisition. Typically, the same hard- 
ware that is used in signal tracking is reconfigured to affect the 
search. During a cold-start, which is a lack of any prior (a 
priori) information about visible GPS signals, acquisition by 
serial search can take upwards of 20 minutes for a very strong 
signal. To acquire weak signals, more data must be examined. 
Approximately 10 times the data record length is required for 
a 10 times weaker signal. Using serial search methods acqui- 
sition times grow quadratically, which results in a 20 -minute 
search increasing to approximately 33 hours, which practi- 
cally means that the signal cannot be acquired at all. Accord- 


5 SUMMARY 

Disclosed herein is a global positioning system (GPS) 
receiver comprising a GPS signal acquisition component 
including a frequency domain correlation module, wherein 
10 said GPS signal acquisition component is adapted to acquire 
said GPS signal with an Effective Isotropic Received Power 
(EIRP) of at least about -180 dBW, and wherein said GPS 
receiver is adapted to operate in a space environment. The 
GPS receiver may further comprise a plurality of GPS signal 
1 5 tracking correlators adapted to track said GPS signals with an 
EIRP of at least about -175 dBW. The GPS signal tracking 
correlators may comprise a Field-Programmable Gate Array 
(FPGA). The GPS signal acquisition component also may 
comprise an FPGA. The GPS signal acquisition component 
20 may be adapted to acquire a GPS signal within about one 
minute for GPS signals with an EIRP greater than about - 1 80 
dBW. The GPS signal acquisition component may be adapted 
to acquire a GPS signal within about one second for GPS 
signals with an EIRP greater than about -160 dBW. The GPS 
25 signal acquisition component alternately may be adapted to 
acquire a GPS signal within about one minute for GPS signals 
with an EIRP greater than about -180 dBW, and to acquire a 
GPS signal within about one second for GPS signals with an 
EIRP greater than about -160 dBW. The GPS signal acqui- 
30 sition component may be adapted to perform a Fast Fourier 
Transform (FFT) correlation process on said GPS signals. 
The GPS signal acquisition component may be adapted to 
acquire a GPS signal without a priori data. 

Also disclosed heroin is a GPS receiver comprising a GPS 
35 signal acquisition component including a frequency domain 
correlation module, wherein said GPS signal acquisition 
component is adapted to acquire a GPS signal within about 
one second for GPS signals with an EIRP greater than about 
-160 dBW, and wherein said GPS receiver is adapted to 
40 operate in a space environment. The GPS receiver may further 
comprise a plurality of GPS signal tracking correlators 
adapted to track said GPS signals with an EIRP of at least 
about -175 dBW. The GPS signal tracking correlators may 
also comprise an FPGA. The GPS signal acquisition compo- 
45 nent optionally may comprise an FPGA. Also, the GPS signal 
acquisition component may comprise a frequency domain 
correlation module. The GPS signal acquisition component 
may be adapted to acquire a GPS signal with an EIRP greater 
than about -180 dBW. The GPS signal acquisition compo- 
50 nent may be adapted to acquire a GPS signal within about one 
minute for GPS signals with an EIRP greater than about - 1 80 
dBW. The GPS signal acquisition component may be adapted 
to perform an FFT correlation process on said GPS signals. 
Alternately, the GPS signal acquisition component may be 
55 adapted to acquire a GPS signal without a priori data for GPS 
signals with an EIRP greater than about -180 dBW. 

Also disclosed herein is a GPS receiver comprising a GPS 
signal acquisition component including a frequency domain 
correlation module, wherein said GPS signal acquisition 
60 component is adapted to acquire said GPS signals with an 
EIRP of at least about -180 dBW without a priori data, and 
wherein the acquired GPS signal is acquired substantially in 
real time. The GPS receiver may further comprise a plurality 
of GPS signal tracking correlators adapted to track said GPS 
65 signals with an EIRP of at least about -175 dBW. The GPS 
signal tracking correlators also may comprise an FPGA. The 
GPS signal acquisition component may comprise an FPGA. 
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The GPS signal acquisition component may be adapted to 
acquire a GPS signal within about one minute for GPS signals 
with an EIRP greater than about - 1 80 dB W. The GPS signal 
acquisition component may be adapted to acquire a GPS 
signal within about one second for GPS signals with an EIRP 
greater than about -160 dBW. The GPS signal acquisition 
component may be adapted to acquire a GPS signal within 
about one minute for GPS signals with an EIRP greater than 
about -180 dBW and to acquire a GPS signal within about 
one second for GPS signals with an EIRP greater than about 
- 1 60 dBW. The GPS signal acquisition component may be 
adapted to perform an FFT correlation process on said GPS 
signals. The GPS signal acquisition component may be 
adapted to acquire a GPS signal without a priori data. 

Also disclosed herein is a method of acquiring and tracking 
GPS signals, comprising the steps of: receiving GPS signals; 
performing a frequency domain correlation process on said 
GPS signals; comparing the correlated GPS signals with a 
predetermined threshold; and acquiring a GPS signal if a 
correlated GPS signal exceeds the predetermined threshold. 
The step of performing a frequency domain correlation pro- 
cess may comprise performing an FFT correlation process for 
a plurality of frequency increments of the received GPS sig- 
nals. The step of acquiring a GPS signal may comprise acquir- 
ing the GPS signal if a correlated GPS signal exceeds the 
predetermined threshold for any of the frequency increments. 
The method of acquiring and tracking GPS signals further 
may comprise receiving GPS signals for a plurality of time 
increments, and accumulating the correlated GPS signals 
over a plurality of time increments to permit the acquisition of 
GPS signals with substantially lower received EIRP. Each of 
the plurality of time increments preferably is about 1 ms. The 
method further may comprise acquiring said GPS signals at 
an EIRP of at least about -180 dBW. The method further may 
comprise tracking said GPS signals at an EIRP of at least 
about -175 dBW. The method further may comprise acquir- 
ing a GPS signal within about one minute for GPS signals 
with an EIRP greater than about -180 dBW. The method 
further may comprise acquiring a GPS signal within about 
one second for GPS signals with an EIRP greater than about 
-160 dBW The method further may comprise acquiring a 
GPS signal with an EIRP greater than about -180 dBW 
within about one minute and acquiring a GPS signal with an 
EIRP greater than about -160 dBW within about one second. 
The method further may comprise acquiring a GPS signal 
with an EIRP greater than about -180 dBW without a priori 
data. 

Also disclosed herein is a method of acquiring and tracking 
GPS signals, comprising the steps: receiving GPS signals for 
a plurality of time increments; performing a frequency 
domain correlation process on said GPS signals for each of 
the plurality of time increments; comparing the correlated 
GPS signals with a predetermined threshold; and acquiring a 
GPS signal if a correlated GPS signal exceeds the predeter- 
mined threshold for any of the plurality of time increments. 
Each of the plurality of time increments preferably is about 1 
ms. The step of performing a frequency domain correlation 
process may include performing an FFT correlation process 
for a plurality of frequency increments of the received GPS 
signals for each of the plurality of time increments. The step 
of acquiring The GPS signal may include acquiring a GPS 
signal if a correlated GPS signal exceeds the predetermined 
threshold for any of the frequency increments for any of the 
time increments . The method further may comprise acquiring 
a GPS signal with an EIRP of at least about -180 dBW. The 
method further may comprise tracking a GPS signal with an 
EIRP of at least about -175 dBW. The method further may 
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comprise acquiring a GPS signal with an EIRP greater than 
about -180 dBW within about one minute. The method fur- 
ther may comprise acquiring a GPS signal with an EIRP 
greater than about -160 dBW within about one second. The 
5 method further may comprise acquiring a GPS signal within 
about one minute for GPS signals with an EIRP greater than 
about -180 dBW and acquiring a GPS signal within about one 
second for GPS signals with an EIRP greater than about -160 
dBW. The method further may comprises acquiring a GPS 
to signal with an EIRP greater than about - 1 80 dBW without a 
priori data. 

Also disclosed herein is a computer program embodied in 
a computer readable medium for acquiring and tracking GPS 
signals, which comprises a receiver module for receiving 
1 5 GPS signals, a correlation module for performing a frequency 
domain correlation process on said GPS signals, a compari- 
son module for comparing the correlated GPS signals with a 
predetermined threshold, and an acquisition module for 
acquiring a GPS signal if a correlated GPS signal exceeds the 
20 predetermined threshold to allow said GPS signal to be 
acquired and tracked. The correlation module may be adapted 
to perform the frequency domain correlation process by con- 
trolling a hardware acquisition component. The hardware 
acquisition component may comprise an FPGA. The com- 
25 puter readable medium may be radiation hardened. The 
receiver module may be adapted to receive GPS signals for 
one or more time increments. The frequency domain corre- 
lation process may comprise an FFT correlation process for a 
plurality of frequency increments of the received GPS signals 
30 for each of the one or more time increments. The acquisition 
module may be adapted to acquire a GPS signal if the corre- 
lated GPS signal exceeds the predetermined threshold for any 
of the frequency increments for any of the one or more time 
increments. The acquisition module may be adapted to 
35 acquire GPS signals with an EIRP of at least about -180 
dBW. The computer program further may comprises a track- 
ing module for tracking GPS signals with an EIRP of at least 
about -175 dBW. The acquisition module may be capable of 
acquiring a GPS signal with an EIRP greater than about - 1 80 
40 dBW within about one minute. The acquisition module may 
be capable of acquiring a GPS signal with an EIRP greater 
than about -160 dBW within about one second. The receiver 
module may receive GPS signals via a hardware receiver. 

The present invention overcomes the problems and disad- 
45 vantages associated with conventional methods and systems, 
and provides improved systems and methods for fast signal 
acquisition and weak signal tracking. 

Other embodiments and advantages of the invention are set 
forth in part in the description that follows, and in part, will be 
50 obvious from this description, or may be learned from the 
practice of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

55 The embodiments of the invention will be better under- 
stood from the following detailed description with reference 
to the drawings, in which: 

FIG. 1 is a flow diagram illustrating a preferred method 
according to an embodiment of the invention; 

60 FIG. 2 is a computer architecture diagram according to an 
embodiment of the invention; 

FIG. 3 illustrates graphical representations of ambiguity 
function cross-sections according to an embodiment of the 
invention; 

65 FIGS. 4 and 5 illustrate graphical representations of per- 
formance curves according to an embodiment of the inven- 
tion; 
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FIG. 6 illustrates graphical representations of receiver 
operating characteristic (ROC) curves according to an 
embodiment of the invention; 

FIG. 7 illustrates a system block diagram according to an 
embodiment of the invention; 

FIG. 8 illustrates a graphical representation of an output of 
a simulation testing an embodiment of the invention; 

FIG. 9 illustrates a dual mode acquisition block diagram 
according to an embodiment of the invention; and 

FIG. 10 illustrates a software/hardware interface schematic 
diagram according to an embodiment of the invention. 

DETAILED DESCRIPTION 

The embodiments of the invention and the various features 
and advantageous details thereof are explained more hilly 
with reference to the non-limiting embodiments that are illus- 
trated in the accompanying drawings and detailed in the fol- 
lowing description. It should be noted that the features illus- 
trated in the drawings are not necessarily drawn to scale. 
Descriptions of well-known components and processing 
techniques are omitted so as to not unnecessarily obscure the 
embodiments of the invention. The examples used herein are 
intended merely to facilitate an understanding of ways in 
which the embodiments of the invention may be practiced and 
to further enable those of skill in the art to practice the 
embodiments of the invention. Accordingly, the examples 
should not be construed as limiting the scope of the embodi- 
ments of the invention. 

As mentioned, there remains a need for a new, fully space- 
qualified weak signal GPS receiver targeted for high altitude 
applications. The embodiments of the invention achieve this 
by providing a space-borne GPS receiver that can operate 
effectively in the full range of Earth orbiting missions from 
LEO to GEO and beyond. Specifically, the embodiments of 
the invention provide a fully space flight qualified GPS 
receiver optimized for fast signal acquisition and weak signal 
tracking. The fast acquisition capabilities provide exceptional 
time to fix performance with no a priori receiver state or GPS 
almanac information, even in the presence of high Doppler 
shifts present in LEO (or near perigee in highly eccentric 
orbits). The fast acquisition capability also makes it feasible 
to implement extended correlation intervals and therefore 
significantly reduce the acquisition threshold of the embodi- 
ments of the invention. This greatly improves GPS observ- 
ability when the receiver is above the GPS constellation (and 
satellites must be tracked from the opposite side of the Earth) 
by providing at least about 1 5 dB of increased acquisition 
sensitivity. One of skill in the art will recognize that the 
present invention, while particularly advantageous for space- 
based applications, from LEO to GEO and beyond, is not 
limited to the particular embodiments enclosed herein and 
can also be applied to other applications known to those of 
skill in the art where weak signal acquisition and tracking or 
fast signal acquisition are desireable. Some exemplary 
embodiments of the present invention are disclosed in the 
paper Luke Winternitz, et al, “Navigator GPS Receiver for 
Fast Acquisition and Weak Signal Tracking Space Applica- 
tions,” Proceedings of the Institute of Navigation GNSS 2004 
Conference, Long Beach, Calif. September 2004, which is 
herein incorporated by reference in its entirety. 

Referring now to the drawings, and more particularly to 
FIGS. 1 through 10, where similar reference characters 
denote corresponding features consistently throughout the 
figures, there are shown preferred embodiments of the inven- 
tion. 
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FIG. 1 is a flow diagram illustrating a method of acquiring 
and tracking GPS signals according to an embodiment of the 
invention, wherein the method comprises receiving (51) GPS 
signals; subdividing (53) the GPS signals into a plurality of 
5 frequency increments; performing (55) a frequency domain 
correlation process on a plurality of the frequency increments 
of said GPS signals; comparing (57) the correlated GPS sig- 
nals with a predetermined threshold; and acquiring (59) a 
GPS signal if a correlated GPS signal exceeds the predeter- 
l o mined threshold for any of the frequency increments of a GPS 
signal. The method optionally further includes tracking (61) 
the acquired GPS signal The step of performing a frequency 
domain correlation process preferably comprises performing 
a Fast Fourier Transform (FFT) correlation process to corre- 
15 late the GPS signals. 

Furthermore, the method may further comprise acquiring 
GPS signals at an EIRP of at least about -180 dBW; and 
tracking the GPS signals at an EIRP of at least about -180 
dBW. Moreover, the method may further comprise fast acqui- 
20 sition of a GPS signal within about one second for GPS 
signals with an EIRP greater than about - 1 60 dBW. In yet 
another embodiment, the method comprises fast acquisition 
of weak GPS with an EIRP greater than about -180 dBW 
within about one minute. One of skill in the art will recognize 
25 that a GPS receiver capable of fast acquisition of both strong 
and weak signals falls within the scope of the present inven- 
tion. In yet another embodiment, the method may further 
comprise acquiring a GPS signal with an EIRP greater than 
about -180 dBW. without a priori data. One of skill in the art 
30 will recognize that the present invention permits acquisition 
of GPS signals more quickly for signals with higher received 
EIRP with no prior knowledge of the GPS signal. 

Additionally, the method may further comprise acquiring a 
GPS signal within about one minute for GPS signals with an 
35 EIRP greater than about -180 dBW. The method optionally 
may further comprise combining the shifted GPS signals 
output by the FFT correlation process with a pre-computed, 
conjugated Discrete Fourier Transform (DFT) output. Addi- 
tionally, the method preferably comprises performing an 
40 Inverse Fast Fourier Transform of the combined GPS signals 
and DFT output. 

The embodiments of the invention can take the form 
including both hardware and software elements. The software 
embodiment includes but is not limited to firmware, resident 
45 software, microcode, etc. Furthermore, the embodiments of 
the invention can take the form of a computer program 
embodied in a computer-readable medium providing pro- 
gram code for use by or in connection with a computer or any 
instruction execution system. For the purposes of this 
50 description, a computer-usable or computer readable medium 
can be any apparatus that can comprise, store, communicate, 
propagate, or transport the program for use by or in connec- 
tion with the instruction execution system, apparatus, or 
device. 

55 The medium can be an electronic, magnetic, optical, elec- 
tromagnetic, infrared, or semiconductor system (or apparatus 
or device) or a propagation medium. Examples of a computer- 
readable medium include a semiconductor or solid-state 
memory, magnetic tape, a removable computer diskette, a 
60 random access memory (RAM), a read-only memory (ROM) , 
a rigid magnetic disk and an optical disk. Current examples of 
optical disks include compact disk — read only memory (CD- 
ROM), compact disk — read/write (CD-R/W) and DVD. One 
of skill in the art will recognize that the foregoing types of 
65 computer readable media are merely exemplary and that other 
types of computer readable media are contemplated by the 
present invention. 
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A data processing system suitable for storing and/or 
executing program code preferably will include at least one 
processor coupled directly or indirectly to memory elements 
through a system bus. The memory elements can include local 
memory employed during actual execution of the program 
code, bulk storage, and cache memories which provide tem- 
porary storage of at least some program code in order to 
reduce the number of times code must be retrieved from bulk 
storage during execution. 

Input/output (I/O) devices (including but not limited to 
keyboards, displays, pointing devices, etc.) can be coupled to 
the system either directly or through intervening I/O control- 
lers. Network adapters may also be coupled to the system to 
enable the data processing system to become coupled to other 
data processing systems or remote printers or storage devices 
through intervening private or public networks. Modems, 
cable modem and Ethernet cards are just a few of the currently 
available types of network adapters. 

A representative hardware environment for practicing the 
embodiments of the invention is depicted in FIG. 2 . This 
schematic drawing illustrates a hardware configuration of an 
information handling/computer system in accordance with 
the embodiments of the invention. The system comprises at 
least one processor or central processing unit (CPU) 10. The 
CPUs 10 are interconnected via system bus 12 to various 
devices such as a random access memory (RAM) 14 , read- 
only memory (ROM) 16, and an input/output (I/O) adapter 
18 . The I/O adapter 18 can connect to peripheral devices, such 
as disk units 11 and tape drives 13 , or other program storage 
devices that are readable by the system. The system can read 
the computer program on the computer readable media and 
execute the computer program to implement the methodol- 
ogy of the embodiments of the invention. The system further 
may include a user interface adapter 19 that connects a key- 
board 15 , mouse 17 , speaker 24 , microphone 22 , and/or other 
user interface devices such as a touch screen device (not 
shown) to the bus 12 to gather user input. Additionally, a 
communication adapter 20 may comiect the bus 12 to a data 
processing network 25 , and a display adapter 21 connects the 
bus 12 to a display device 23 which may be embodied as an 
output device such as a monitor, printer, or transmitter, for 
example. 

Another embodiment of the invention provides aGPS 
receiver adapted for use in a space environment. In yet 
another embodiment, the GPS receiver may be adapted for 
operation in LEO orbit, a GEO orbit, or a variety of orbital 
conditions above or below GEO orbit. The GPS receiver may 
be fully space flight qualified using radiation-hard parts. 
Another embodiment preferably implements methodologies 
or algorithms in field-programmable gate arrays (FPGAs) to 
allow for easy modification, upgrading, and customization. 
Alternately, the method of the present invention may be 
implemented in one or more ASICs. One of skill in the art will 
recognize that a variety of hardware combinations fall within 
the scope of the present invention. 

Another embodiment provides a GPS receiver which can 
acquire and track signals at an EIRP of -180 dBW. Another 
embodiment of the invention provides a GPS receiver capable 
of acquiring GPS signals within about one second for strong 
signals (with an EIRP greater than -160 dBW) and one 
minute for weak signals (with an EIRP greater than -180 
dBW). One of skill in the art will recognize that the present 
invention permits acquisition of GPS signals more quickly for 
signals with higher received EIRP with or without a priori 
knowledge of the GPS signals. 

The embodiments of the present invention provide perfor- 
mance improvements over the conventional receivers: ( 1 ) 
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signals can be reliably acquired and tracked down to a 
received EIRP of about -180 dBW, which is at least a 15 dB 
improvement in sensitivity; and (2) Fast Fourier Transform 
(FFT) based acquisition algorithms provide extremely short 
5 acquisition times even for the weakest signals. 

In order to make use of a GPS signal coming from a 
particular satellite, a GPS receiver first acquires and then 
tracks the signal. Again, the acquisition process is a much 
more computationally demanding task than the tracking pro- 
to cess, requiring a search across a three-dimensional space of 
unknown time delay, Doppler shift, and satellite number. It is 
also a limiting factor for sensitivity for the common GPS 
receiver architecture. The acquisition methodologies imple- 
mented in the embodiments of the invention may be adapted 
1 5 from the weak GPS signal acquisition technique provided by 
Psiaki, M. L., “Block Acquisition of Weak GPS Signals in a 
Software Receiver,” ION GPS 2001, 11-14 Sep. 2001, Salt 
Lake City, Utah, the complete disclosure of which, is herein 
incorporated by reference in its entirety. Psiaki presents a 
20 detection statistic and hypothesis testing scheme that enables 
the acquisition of very weak GPS signals with no a priori 
knowledge; in particular, no knowledge of the GPS data mes- 
sage bits. An embodiment of the present invention employs 
the methodologies taught in Psiaki using specialized hard- 
25 ware implemented in radiation-hardened FPGAs. In another 
embodiment, the acquisition process is capable of running 
substantially in real-time. For the purposes of this embodi- 
ment, substantially in real time means that the signal param- 
eters of the acquired signal are valid for tracking purposes 
30 when acquired. In another embodiment, the GPS receiver is 
capable of achieving the desired sensitivity with acquisition 
times of only a few seconds. 

The received, down converted and sampled GPS signal can 
be preferably modeled as: 

35 

y k =s k +n k k = 1, 2, ... (1) 

s k = — “W~7^ — -\e^ kTs = Ad k c k e^ kTs 

V 1 + re / V 1 +re / 

40 

= \A\e> 9 = cr / T s eP , n k ~ CN( 0, v 2 ) 

V No 


45 Computationally, the received signal is y k , comprising the 
sum of the GPS signal, s k , and complex white Gaussian noise, 
% (real and imaginary part of % are N(0,cr 2 /2) and indepen- 
dent). Here, d(.) is the 50 bps GPS data message and c(.) is the 
1.023 Mcps, 1 ms (1023 chip) periodic Pseudo Random 
50 Noise (PRN) spreading code, both taking values in {-1 ,+l }. 
Moreover, A comprises the complex GPS signal amplitude/ 
phase, which is related to the received C/N 0 as shown with T^, 
which comprises the sampling period in seconds. Moreover, 
co connotes the Doppler shift in rad/ s and r\ =to/ 1.57542xl0 9 is 
55 the unitless Doppler compression factor on the PRN code. 
Additionally, t represents the time delay in the seconds mod- 
ule and is approximately 0.00 1 sec. Furthermore, c k and d k are 
the appropriately compressed and delayed code and data 
samples (c^(t) is denoted when it is desired to make depen- 
60 dence on t explicit). The noise can be modeled as white 
because it is assumed that the signal is bandpass filtered to 
capture only the mainlobe and then sampled at the Nyquist 
rate for the resulting bandwidth. 

For simplicity, it preferably is first assumed that there is 
65 only one known satellite signal potentially present and the 
data modulation d(.) is ignored. Then, GPS signal acquisition 
comprises determining whether this signal is present and, if 
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so, at which point, (co,x). Given a finite set of samples, y=(y 0 , is declared at that particular (oo,x). In the serial search, new 

y ls . . . ,y N _i) (s, n defined analogously), it is desired to test data is used at each grid point. This amounts to performing an 

between the following hypotheses, and if H 1 is accepted, independent binary test at each grid point with (oo,x) known 
(oo,x) is preferably estimated as: (thus no need to estimate them). In fact, if phase is treated as 


H o :y=n < ^ > \A\=0 

H l :y=s+n <f> \A\>0 (2) 

The signal, s, may depend on a number of unknown param- 
eters (A,oo,x) which lie in some known space. The true values 
of these parameters are denoted as: (A.),co.,x.). Because of 
these unknowns, no clear optimal test generally exists under 
any common criteria. One approach is to use a Generalized 
Likelihood Ratio Test (GLRT). This procedure will usually 
produce the optimal test if it exists and otherwise generally 
gives effective sub-optimal testing procedures. The GLRT 
prescribes the following test: 


„ rD sup {Ar(0rT )My,A,(o, T ) (3) 

GLR = > y 

fo(y) 

Here, f z (.) denotes the joint probability density of the data 
under hypothesis H : . That is, the ratio of the corresponding 
joint probability densities of the data under each hypothesis is 
computed, with the observed data substituted in (the likeli- 
hoods). Then, is decided if this ratio exceeds some thresh- 

old, which is preferably chosen to fix or bound the false alarm 
probability, using the Ney man- Pearson (NP) criterion. 

Unknown parameters preferably are replaced by their maxi- 
mum likelihood (ML) estimates under The ML estimate 
of A may easily be determined (as a function of to andx) to be: 

l £} ( 4 ) 

A ML (a), t) = jj2j y k c ke MTs 

k=0 

\A\ml = \Aml\, 0ml = LAml 


The maximization of Equation (3) over (co,x) may be per- 
formed by a grid search over the range of allowable values: 
(oo,x) in [oo mzz7 , co mz2X]x[ 0, 0.001]. It can be shown by substitut- 
ing in the corresponding multivariate Gaussian densities, that 
Equation (3) reduces to the familiar test on the correlation 
statistic (scaling term included to normalize the noise): 


max^r) 13^)001 > y (5) 

ZferlOO = 

' k=0 

That is, one looks for the maximum magnitude- squared, 
over (oo,x), of the correlation between the input and a unit 
amplitude, zero -phase replica of the signal and compares it 
against a threshold. This is equivalent to comparing the full 
ML estimate of the signal amplitude against a threshold, 
which is intuitively satisfying. 

The test described above preferably utilizes the calculation 
of correlation at each (oo,x) using the same set of input 
samples. This is referred to as the parallel search. In contrast, 
most conventional GPS receivers typically employ a serial 
search. The local signal generator serially steps through the 
(oo,x) grid and computes a correlation with the input samples, 
as the data streams in. When a threshold is crossed, a detection 


5 random and uniformly distributed in (0, 2 ji), this is a preferred 
(Uniformly Most Powerful for I Al >0) test under the NP and 
Bayes criteria. In this case, there is no maximization and the 
statistics of Z, from Equation (5), completely, specify the 
performance of the test. The false alarm probability prefer- 
10 ably is set for the individual tests and thus is actually the false 
alarm rate. This rate can be traded off against the detection 
probability. In the serial search, false alarms slow the acqui- 
sition process, but are not devastating because the tracking 
eventually fails when initialized with a false signal. After the 
15 failure, acquisition can continue where it is left off. 

In the parallel search, a single block of data is used to 
maximize IZI over the test grid. Although the same computa- 
tions are performed as in the serial case (only on a fixed block 
of data rather than new blocks), the problem is statistically 
20 very different. In this case, the test statistic is maxi Z I whose 
exact distribution is very difficult to obtain since there is a 
dependence between the correlations across the search grid. 
The “(tt),x) known” statistics, which hilly characterize the 
serial test, can be used in the parallel case to obtain bounds on 
25 the false alarm probability as well as the approximate perfor- 
mance. 

One method for controlling in the parallel search, is 
via the union bound. Here, one preferably divides the desired 
overall false alarm probability by the number of grid points 
30 and uses that to set the threshold assuming (oo,x) known. In 
another embodiment, one can assume: 


Pr(max|Z(k>, r)| > y) % Pr(max|Z(k/, r')| > y) (6) 


where the (co^x’) are the mutually independent points on the 
grid. There are approximately 

40 K We/z =1023*DopplerRange*NT, of these independent 
points. With independence, one can solve for the exact thresh- 
old for a given P^. Using P^/K indep in the (co,x) known test 
provides a very good approximation to the desired threshold. 

False alarms are generally less tolerable in the parallel 
45 case. There is a range of threshold values where there i s a very 
high probability that at the true signal point, Z(oo 0 ,x 0 ) will 
exceed the threshold, and so the signal should eventually be 
detected by the serial search. However, there is also a rela- 
tively large false alarm rate and a large probability that Z(oo 0 , 
50 x 0 ) is not the largest correlation over the grid. One way to look 
at it is that the correct decision, H l5 is made, but the estimation 
of (oo,x) is bad. This is referred to as a “Type III” error, as 
opposed to false alarms under H 0 or “Type I” errors, and 
missed detections under H 1 or “Type II” errors. Another 
55 approach to this problem is as an M-ary test with H (cox) : 
“signal present at (co,x)”, wherein similar results can be 
obtained. 

One of the problems with a serial search is that the serial 
search process, where new data is used in each correlation, 
60 takes too long for weak signals. For example, using N corre- 
sponding to 1 ms typically allows for a best-case acquisition 
of signals of about 35 dB-Hz, which corresponds to an EIRP 
of approximately -165 dB W, (with many false alarms). In this 
case, the time to first fix performance (TTFF) could be on the 
65 order of about 30 min. If one desires a 1 0 times improvement 
in sensitivity, one could increase the N to 10 ms, but in 
addition one also has to increase the fineness of the grid in the 



US 7,548,199 B2 


11 


12 


frequency dimension by another factor of 10. This implies a 
new TIFF on the order of about 3,000 min. Accordingly, the 
search preferably is parallelized. In either case, a rather large 
number of correlations need to be quickly computed. 

The following results from computing Equation (4): 5 


Z(to,T)(y) 




2 x 

j ^~2 A^c k (T 0 )c k )k T - + n= 


■M 


NT s ej( d+ ») + n= 



N 0 


|/?(Ar, Au))\ 2 


n= ~ CN( 0, 2) 


(7) 


10 


15 


where (Aoo, Ax)=(oo-co 0 , x-x 0 ) comprise the errors in estimat- 
ing (oo,t), and 0 (tilde) comprises the phase of R(Aco,At). This 2 o 
R(Aoo,At) (which approximately depends only on Ax) is 
sometimes referred to as the ambiguity function, and it gen- 
erally specifies the necessary fineness of the grid upon which 
the maximization of Equation (4) is performed. Uncertainty 
or error in the estimates of to and x result in a reduced mean of 2 5 
the correlation statistic and can be viewed as a decrease of the 
effective input C/N 0 . The cross sections of R(.,.), which are 
plotted in FIG. 3 are given by the following equations: 


divided down by 1023 (appropriate for search across all code 
delays and a 1 kHz Doppler range) are shown in FIG. 4 , with 
P^, vs. coherent integration time (oo,x) known and P^=0.05/ 
1023. Here, a factor of 10 increase in integration time gener- 
ally results in a 10 dB increase in sensitivity. 

Unfortunately, the 20 ms GPS databit limits how long the 
coherent integration can run (other effects such as signal and 
clock dynamics limit this time as well but on a longer scale 
than 20 ms). Correlating over a databit transition may result in 
unpredictable cancellation of correlation power. Even using 
20 ms generally requires knowledge of the bit start times. The 
bit start time could be added to the problem as an unknown 
parameter and one could further maximize Equation (3) over 
the additional dimension. Generally, a single 20 ms coherent 
correlation generally will not reliably detect the desired 25 
dB-Hz signal (corresponding to an EIRP of approximately 
-175 dBW), accordingly, what is needed is a way to extend 
the correlations beyond the databits. 

One embodiment capable of solving the problem is the 
“half-bits” method, which collects consecutive 10 ms corre- 
lations. If one looks at this collection as two alternating sets, 
then, one is guaranteed to be free of transitions. A collection 
of M 10 ms (or any arbitrary L ms not necessarily equal to 1 0) 
coherent correlations that individually avoid databit transi- 
tions, {Z ((OT) m5 m=l,2 . . . L}, can be combined “non-coher- 
ently” by taking the magnitude squared and then summing. 
This removes the problematic databit. 


|/?(Ar, 0)| 2 » 



sm(Au>T s N / 2) 2 
Nsm(Aa)T s / 2) 


30 

( 8 ) 


M-l 
m= 0 

Q* = maX( W)r) Q {oj , t) 


( 10 ) 


It is desirable to keep I Acol <1/(4 NT S ),\ restricting loss along 
a frequency axis to approximately 0.2 dB and to keep IAI<14 
chip, which restricts the loss to approximately 2.5 dB along 
the code-delay axis. In FIG. 3, Ao)=0 in the upper graph and 
Ax=0 in the lower graph. 

Assuming the (oo,x) are known or fixed, and the GPS data 
message is constant over the block of input samples of inter- 
est, then the statistics of Z are readily determined to be: 



z^rtOO = yfte* +rh ~ cn(4x^, 2 ) 

Ho ■■ IZ^riOOl 2 ~ X|( 0) 

Ml : \Aoj,r)(yf ~ x 2 2 (A) 

Under H 0 the test statistic is a chi-squared random variable 
with two degrees of freedom and under H 1 the test statistic is 
a second degree non-central chi-squared with non-centrality 
parameter X, which is a product of twice the effective C/N 0 
and the integration time, NTs. 

These statistics fully characterize the serial search and can 
be used to bound the false alarm probability for the parallel 
search. They also provide the probability that the true signal 
will cross the threshold under H l5 but do not account for the 
“Type III” event that there is a larger correlation at another 
( 00 , x) coordinate. When the overall false alarm probability is 
kept low, however, the “Type III” event is rare and the curves 
shown in FIG. 4 provide an approximate performance. Spe- 
cifically, performance curves for the parallel search with 


Q* is then used as the test statistic. This is the “P lan ” detec- 
tion statistic described by Psiaki. This seems reasonable, and 
can also be arrived at as the Locally Most Powerful (LMP) 
test under certain assumptions. This combined L ms coherent 
and M block non-coherent integration is referred to as the 
“L/M integration”. 

Another reason to keep the coherent integration period 
short is because of the structure of the ambiguity junction 
along the frequency axis. Correlations preferably can be com- 
bined non-coherently for as long as desired without the need 
for increasing frequency resolution. For 1 0 ms coherent inte- 
gration, a 25 Hz grid spacing will generally result in a worse 
case loss of only approximately 0.2 dB. 

Again, for a given (co,x) and assuming the databit transi- 
tions have been avoided, the statistics of Q may be determined 
to be: 

H o :Q(^)(y)~X 2M 2 (0) (11) 

Performance curves for the parallel search with P FA divided 
down by 1023 are shown in FIG. 5, with P^ vs. total integra- 
tion time for 10M integrations with (00 ,x) known and 
.05/1023. The gains from non-coherent integration 
come more slowly than from coherent integration. Sensitivity 
increases roughly with the square root of integration time. A 
10 times longer integration results in only approximately a 5 
dB increase in sensitivity. However, a 25 dB-Hz signal can be 
very reliably acquired (>0.9 probability) by summing the 
squared magnitudes of only 1 0 ms coherent correlations. This 
is only one tenth of a second of data. Furthermore, 20 dB-Hz 
and lower are reachable with not unreasonably long data 
records. 


35 


40 


45 


50 


55 


60 
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As a check on the reliability of these results, FIG. 6 shows 
the results of a small simulation study which examines the 
Receiver Operating Characteristic (ROC, vs. P^) for a 25 
dB-Hz GPS signal with zero Doppler assumed. In FIG. 6 the 
ROC curves are shown for ROC for 10/10 integration, 25 5 
dB-Hz signal, with co=0 known, t unknown, P FA divided 
down by 1023. The observed curve in FIG. 6 shows good 
agreement with the theoretical result. The Type III error 
deviations are also clearly seen to be a real phenomenon but as 
mentioned, they are rare when the p« is low and thus the 
threshold is set high. 

One of the primary effects of interfering signals is that the 
cross-correlation peaks cause the number of false alarms to 
increase for a given threshold level. One solution is to 15 
increase the thresholds to control the P FA and correspondingly 
increase the integration intervals to recover P^. This will 
work when the power difference is moderate, but may fail for 
disparities greater than about 20 dB, at which point the acqui- 
sition of the weak signal is difficult without somehow dealing 20 
with the interfering signal. One possible solution is through 
careful selection among the available antennas on the space- 
craft or through the use of antenna phasing to help reject the 
unwanted strong signal jammer. Alternatively, one can try to 
acquire and track the strong signals first and then cancel them 25 
out of the input before attempting the weak signals. During 
tracking, precise estimates of (cl),t,0 ) are maintained and all 
that is additionally needed is the signal amplitude, I Al , whose 
ML estimate is given in Equation (4). It is possible to use this 30 
in online processing and may be incorporated in the embodi- 
ments of the invention. 

In one embodiment, the receiver can quickly and reliably 
acquire and track signals down to received EIRP of approxi- 
mately -180 dBW and lower by employing special hardware 35 
that computes the Q-statistic described above. In another 
embodiment, the receiver preferably operates autonomously. 
That is, it does not require external data aiding or any other a 
priori information such as a current estimate of time, a recent 
GPS almanac or a converged navigation filter estimate of the 40 
receiver dynamics, etc. If such a priori information is avail- 
able, then the embodiments of the invention may be able to 
make use of it to reduce acquisition time, but it is not assumed 
to be available and is not required. 45 

In yet another embodiment, the GPS receiver of the present 
invention preferably operates substantially in real-time. Data 
is double buffered up front in 1 ms blocks and processed as it 
comes in. One could use large up front storage to buffer the 
entire record of sampled data needed for acquisition. Process- 50 
ing can then occur with relaxed elite constraints, only needing 
to be sufficiently fast such that the (oo,t) estimate will still be 
valid when processing finishes. In yet another embodiment, 
the receiver operates in a faster than real-time mode to catch 
up with the streaming input, or simply begins tracking on the 55 
most recent samples. 

Yet another embodiment uses Discrete Fourier Transform 
(DFt) based computation of the 1 ms correlations. When the 
FFT methodology is used to compute these DFTs, large com- 6f) 
putational savings are achieved. Specifically, when an 
N-point DFT is computed, the savings are 0(N/log 2 (N)). This 
method preferably calculates all code delay correlations in Vi 
chip increments using a single 1 ms block of data. This 
technique optionally may be used for Direct Sequence Spread 65 
Spectrum (DSSS) signal acquisition. Generally, the DFT of a 
signal x=(x 0 , x 1 . . . x^)'; preferably is given by: 


N-l 2x 

X k =Y J X«e- Jwkn =DFT[x] k 

n = 0 


(12) 


The “circular-correlation property” gives the following 
identity (where ° connotes point by point multiplication of the 
vectors): 


N - 1 


Z n — ) ' XkC^+rijmodN — IDFT\X °C ] w 

k= 0 


(13) 


To apply this to the current problem, a 1 ms block of the 
baseband downconverted input signal is used as the x and the 
DFT of the code sequence as C, which can be computed 
offline, and the following is computed (t now represents an 
integer offset of the code between 0 and N-l): 


ZtariOO = L b'i <r ‘ M7 4^ 


(14) 


N-l 

= X {yte~ MTs ]C(k+T )m *w = IDFT[X °C*] T 

k=0 


This gives the 1 ms correlation. The effect of Doppler shift 
on the code sequence is ignored for now and is further dis- 
cussed below. To get a longer correlation, L consecutive 1 ms 
block correlations are added: 


LN- 1 (15) 

Zinviy) = Yj y^ c ^ e ~ MTs 

k=0 

L - 1 (N+l)/-l 

= ) j ) j \.yk& j S l c (£+r)mo<iV 
/=0 k=Nl 

L - 1 

= lDFT[X t oC*] T 
1=0 


This calculation can be repeated and the squared magnitudes 
summed to get the desired Q-statistic. In one embodiment, the 
acquisition hardware implements this sequential process. In 
another embodiment, the sequential process is implemented 
in software. 

With this method, the Q-statistic is calculated for the entire 
code dimension based on a single block of input data. Fur- 
thermore, if many FFT s can be performed in bits, then one can 
search across the Doppler dimension as well. To search the 
frequency dimension, this method preferably is repeated for 
each frequency on the search grid at whatever frequency 
granularity desired. In this way, the entire (co,x) grid (at least 
over a range of Doppler frequencies) can be searched in 
parallel. 

The embodiments of the invention may utilize different 
methods to reduce the number of DFT operations required. 
One method preferably uses the well-known frequency shift- 
ing property of the DFT. A reference (coarse) Doppler shift 
preferably is removed from a 1 ms block of the input signal 
and its DFT is computed. The DFT vectors of all multiples of 
1 kHz offset from this reference are then obtained by shifting 
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the reference DFT vector. This technique reduces the number 
of forward DFTs needed, but not the number of Inverse Dis- 
crete Fourier Transforms (IDFTs). Another method explains 
how to interpolate from a coarse grid of frequencies with 
spacing on the order of hundreds of Hz to the needed fine grid 5 
spacing of approximately a factor of 10 finer. This interpola- 
tion onto the fine grid is similar to the “Post-Correlation FFT” 
methods that are becoming popular for GPS acquisition. 
However, a 1 ms coherent dump time preferably is used, as 
opposed to the fractional ms dump conventionally used, and to 
rather than using an FFT, which computes the Discrete Time 
Fourier Transform (DTFT) ordinates at the Fourier frequen- 
cies, l/(LNTj, the embodiments of the invention preferably 
compute a few DTFT values at the desired frequencies. 

The first equality in Equation (14) holds when c(.) is 15 
exactly N-periodic. An appropriately Doppler compressed/ 
expanded version of the code that would not be exactly N-pe- 
riodic prerferably is correlated with the input. Exemplary 
methods to solve this problem include resampling the input in 
each Doppler bin to enforce exactly 1 PRN period in N 20 
samples. In this case the N-periodic c(.) is the correct replica. 
Alternately, one could implement a delay correction applied 
in the frequency domain prior to the IDFT. 

A preferred approach is to initially ignore the issue and 
correct for the ill effects later. One embodiment preferably 25 
employs a re-mapping of the stored correlation grid. The 
effect of the code Doppler is that the apparent correlation 
peak drifts through the correlation grid in the delay dimension 
at a rate determined by the Doppler bin under test. It is 
sufficient to correct this issue only upon accumulation into 30 
memory approximately every 20 ms (via re-mapping back to 
the appropriate code bin). 


Drift = ^ chip J dopp chips/sec 
JLl 


Assuming an 8 kHz maximum Doppler magnitude, appropri- 
ate for GEO, this drift is limited by 0.1039 chips in 20 ms. 40 
This actually serves to average out the worst (and best) case 
loss caused by code misalignment for very long acquisitions. 
Using Vi chip search spacing, rather than 2.5 dB worst-case 
loss, an average loss of approximately 1 .2 dB is achieved for 
long integrations affected by moderate Doppler shifts. 45 

In yet another embodiment, the present invention offers an 
extremely fast strong signal mode. In this mode, the acquisi- 
tion is based on a single 1 ms block of input samples. This 1 
ms block optionally is buffered and then FFT based correla- 
tions preferably sweep across the entire Doppler space (Dop- 50 
pier removal) at the desired granularity. A frequency-shifting 
device preferably is used in this mode, but interpolation onto 
the fine grid is not needed. The maximum correlation value of 
this operation preferably can be determined without the need 
to store the entire grid. Thus, the need for a large amount of 55 
high bandwidth memory preferably is substantially elimi- 
nated in this mode. 

Using the maximum over the correlation grid as the detec- 
tion statistic will necessarily limit the sensitivity of this mode 
to received EIRP of about -160 dBW. The threshold prefer- 60 
ably is set rather high to limit the “Type III” errors. Other 
conventional receivers may be able to acquire signals at about 
-165 dBW based on the same 1 ms correlation but they tend 
to do this at the cost of very high false alarm rates (i.e., 
approximately 0.17%). 65 

Another embodiment of the present invention offers a 
weak signal mode. In this embodiment, 10 ms coherent inte- 


grations are combined (separated by a 10 ms delay) non- 
coherently by summing the squared magnitudes as in the 
definition of Q in Equation (10). This can be performed 
indefinitely, limited only by signal and clock dynamics. In 
practice M=10 will acquire the -175 dBW signal reliably and 
M=100 can achieve -180 dBW. 

The weak signal mode utilizes two separate memories. The 
memory preferably stores the 1 0 individual 1 ms correlations 
for each coarse Doppler and the accumulation of the Q sta- 
tistic. Preferably, two Q-statistics are collected, one for each 
set of alternating 10 ms blocks. However, only one preferably 
is computed at a time. The “off’ 10 ms preferably is used for 
computation. The second set preferably is checked only if the 
first set fails to produce a detection. Since the coherent inte- 
gration period is 1 0 ms, a 25 Hz Doppler spacing is preferably 
used, which restricts C/N 0 losses to approximately 0.2 dB. 
Due to the fine granularity, this mode uses both the frequency 
shifting properly of the DFT and interpolates onto a fine grid 
of frequencies to reduce the number of FFT operations 
needed. 

The maximum Doppler range that can be searched at once 
is limited by the number of FFT operations that can be com- 
pleted in 1 ms and the available bandwidth to the off chip 
memories. This is further described below. To cover the fill 
Doppler range, the frequencies are searched sequentially in 
this maximum Doppler block size. Arbitrary T/M modes are 
available as well, but 1/1 (strong) and 10/M (weak) are deter- 
mined to be the most useful. For example, if data aiding is 
available then L>10 may be desirable. 

In the weak signal mode, the ( 00 , t) estimate provided by the 
acquisition module may not be sufficiently accurate to initial- 
ize tracking, particularly in the 00 dimension. Furthermore, 
determination of the location of the GPS databit transition is 
significant when tracking weak signals. When the acquisition 
terminates with a successful detection, the estimated ( 00 , x) is 
used to initialize a tracking correlator/channel. Open loop 
1-PRN (~1 ms) period accumulations are collected for a 
specified amount of time; e.g. J-PRN periods, during which 
the following statistic preferably is calculated to determine an 
estimate of the bit transition time; i.e., to get bit lock (K=floor 
(J/20)): 


.2 


1 20 * + 19 +/ 

z ^ 

7=20*+/ 


(17) 


This assumes that the initial frequency estimate is sufficiently 
good enough so residual frequency error does not cancel out 
the above 20 ms correlations. In another embodiment, I Aco I is 
preferably verified to be less than 12.5 Hz before attempting 
to compute Equation (17). 

To resolve frequency more finely, one can process the same 
type of {ZJ post-correlation sequence. When the replica 
phase is kept continuous over K consecutive T ms coherent 
integration blocks, the individual correlations can have the 
preferred form: 

Z k =bd^ A,x,k(LNTs) +n k (1 8) 

Here, d k is the GPS databit and b is essentially a complex 
constant independent, of k. This being so, if one were to 
square each Z k (to remove the databit), then one acquires a 
pure complex exponential plus (non-Gaussian) noise, for 
which the maximum of the periodogram is an often used 
estimator of frequency. One can use an FFT engine to com- 
pute the periodogram and maximize it to get an estimate of 
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2Aoo. Experimentally, this has been shown to work using L=1 
ms, for signals even below 25 dB-Hz with J=2048, giving a 
resolution of <1 Hz. for Aoo. For smaller J, one may zero-pad 
the Z-squared sequence. This assumes the frequency rate is 
sufficiently small such that the residual Doppler is roughly 
constant over J*L-ms. Depending on the application, the two 
calculations can be performed simultaneously, or first Equa- 
tion (17) then Equation (18) is performed, or vice-versa. 

Another embodiment employs standard frequency locked 
loop/phase lock loop/delay lock loop (FLL/PLL/DLL) track- 
ing methods using 1 ms correlations for strong signals and 
with correlations extended to 20 ms for tracking weak signals. 
PLL tracking down to approximately 23 dB-Hz is achievable 
with 20 ms correlations and an oscillator. Phase tracking is 
known to be the weak link as compared with carrier frequency 
code tracking. The theoretical bit error probability at 
C/N 0 =25 dB-Hz is approximately 10 -4 for perfect coherent 
demodulation of binary phase shift keying (BPSK) signals, so 
data demodulation can be performed relatively well. Much 
below 25 dB-Hz, using traditional FLL/PLL/DLL tracking, 
reliable data demodulation becomes difficult and phase track- 
ing begins to fail. Moreover, bit error probabilities can be 
improved by averaging many cycles of the repeating GPS 
data message. 

As illustrated in FIG. 7 , the receiver 100 provided by one 
embodiment of the present invention comprises an antenna 
101 operatively connected to an analog front end component 
102, which is operatively connected to a bank of hardware 
tracking correlators 103 controlled by a general purpose 
microprocessor 105 . In addition, the analog front-end com- 
ponent 102 also operatively connects to a specialized acqui- 
sition module 104 that rapidly calculates the long-term detec- 
tion statistic, Q*, and is controlled by a memory device 106 
and the microprocessor 105 . 

All of the GPS specific hardware shown in FIG. 7 prefer- 
ably is implemented in Very High Speed Integrated Circuit 
(VHSIC) Hardware Description Language (VHDL) to target 
radiation-hardened FPGAs. Alternately, an Application Spe- 
cific Integrated Circuit (ASIC) or other hardware implemen- 
tation. In yet another embodiment some or all of the hardware 
components depicted in FIG. 7 may be implemented in soft- 
ware without departing from the scope of the present inven- 
tion. 

Simulated software models of the above methodologies 
were used to evaluate the performance and implement modi- 
fications. Other models were designed to emulate the hard- 
ware and were used as a guide for the hardware design. 
Specifically, sampling rate, analog-to -digital converter 
(ADC) bit depth, and all other data path bit-depths, most 
importantly for the FFT, were determined. For hardware 
implementation, these rates and depths are preferably held to 
be as small as possible without sacrificing too much perfor- 
mance. 

Experimentally, the sampling rate was chosen to be 
approximately 2.048 Msps. This is just over the Nyquist rate 
for capturing the main lobe of the C/A code (using complex 
samples) and gives 2048 samples in the 1 ms nominal code 
period. This allows for the use of the efficient radix-2 FFT 
algorithm. This embodiment yields a Vi chip resolution for 
the acquisition that amounts to a 2.5 dB worst-case loss of the 
effective C/N 0 . This loss may be significant and so 4.096 
Msps is also considered. However, doubling the sampling rate 
also tends to double the memory size and bandwidth require- 
ment, which may be deemed an unacceptable cost. 

For the simulations of certain embodiments of the inven- 
tion, ADC depth was chosen to be 8 -bits in the acquisition 
module in order to maintain a large dynamic range. It is 
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determined that 10-bits is the preferred bit-depth for the FFT. 
Simulations of “acquisition-margin” (ratio of largest to next 
largest correlation) vs. bit-depth show that sensitivity begins 
to fall off with fewer than 1 0-bits, while very little is gained by 
5 having more. FIG. 8 shows the output of one such simulation, 
which illustrates the acquisition margin vs. FFT bit dept for 
25-29 dB-Hz signals, 10/10 integration, with 2,000 runs aver- 
aged. 

In yet another embodiment, the acquisition hardware, with 
to reference to FIG. 9, preferably accepts as input: a satellite 
number, a Doppler range, Doppler search granularity, and 
coherent andnon-coherent integration times. Then, the acqui- 
sition hardware 104 (of FIG. 7 ) preferably calculates the 
Q-statistic using the following sequence of steps, again with 
15 reference to FIG. 9: 

1 . The acquisition FPGA receives input samples from the 
ADC. These are streamed into a 2 ms double buffer 107 . 

2 . Samples from the ready buffer 107 enter the Doppler 
removal block 108 which removes the nominal carrier and 

20 residual (coarse) Doppler under test. 

3 . Next the FFT operation 109 operates on 2048 samples 
and the output is stored in a 2048 -point frequency shift buffer 
110 . 

4. Shifted versions of this output are multiplied 111 against 
25 the pre-computed, conjugated Discrete Fourier Transform 

(DFT) of the desired code 112 and then Inverse Fast Fourier 
Transformed (IFFT’d) 113 to complete the 1 ms correlation 
114 . Each one-sample shift corresponds to a 1 kHz offset 
from the originally removed Doppler. 

30 Strong signal mode acquisition finishes by sweeping 
across the Doppler range by repeating steps 2-4 and then 
determining the maximum magnitude (and coordinates) 115 
on the 1 ms correlation grid (including the peak correlation 
detector 116) which it reports back to the microprocessor 105 . 
35 The weak signal mode finishes in the following way: 

5. Each 1 ms vector preferably is stored to the short-term 
memory 1 17 by repealing steps 3 -4 for each coarse frequency 
bin in the search range. Programmable rounding is applied to 
allow 8-bit representation. 

40 6 . At the end of 10 ms, during the “off’ 10 ms, for each 

coarse frequency, the 1 0 corresponding 1 ms correlation vec- 
tors 118 are read from the short-term memory 117 . The inter- 
polation from the coarse to the fine grid is applied at this point 
to achieve the desired fine frequency resolution. 

45 7 . The squared magnitudes 119 of the correlations 118 are 

computed 120 and the code Doppler correction 122 is applied 
before it is accumulated into the long-term memory 129 . 
Rounding here is used to maintain 16-bit words in the long- 
term memory 129 . One of skill in the art will recognize that it 
50 is possible to combine short-term memory 117 with long- 
term memory 129 into a single memory device or a plurality 
of memory devices without departing from the scope of the 
present invention. 

At the end of the accumulations, the maximum over the 
55 long-term grid, Q*, preferably is determined with the peak 
correlation detector 116 and reported to the microprocessor 
105 for threshold comparison. 

Some potential bottlenecks limit the performance of the 
acquisition module 104 (of FIGS. 7 and 10 ) including: the 
60 bandwidth to off-chip static random access memory (SRAM) 
(for accumulation of the correlation grids) and speed of the 
FFT operation 109 (of FIG. 9 ), which is ultimately deter- 
mined by the size of the FPGA. Running this methodology in 
real-time requires a relatively large bandwidth to the off chip 
65 SRAM (not shown). To provide for this bandwidth, a 64 -bit 
bus (not shown) connects the SRAM to the acquisition FPGA. 
For example, if the 64 -bit bus operates at 66 MHz, this pro- 
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vides a 528 MB/s bandwidth. Weak signal mode requires 
approximately 32.8 MB/s bandwidth per 1 kHz Doppler 
search block and thus can cover about 16 kHz at once if 
limited only by memory bandwidth. 

Another potential bottleneck is running the FFT s with suf- 5 
ficient speed. In one embodiment, the hardware has the capa- 
bility of performing 23 FFTs/ms. Assuming a 250 Hz coarse 
Doppler spacing, this implies a maximum one pass Doppler 
coverage of 5.5 kHz. One embodiment of the present inven- 
tion utilizes 4 parallel butterfly adders (not shown) to imple- 10 
ment the FFT, limited by the expected utilization of the flight 
FPGA. If the usage estimates are proven to be too conserva- 
tive, the number of butterfly adders could be increased, pro- 
viding improved FFT speed. 

The tracking FPGA preferably comprises a standard block 1 5 
of hardware correlators 103 (of FIG. 7). However, to improve 
FPGA usage efficiency, rather than processing samples at the 
sampling rate, data preferably is stored in a first-in-first-out 
(FIFO) manner and processed by time shared hardware run- 
ning at a much higher rate than the sampler. The specific 20 
implementation preferably includes three time-shared corr- 
elator blocks that give 12 channels each, for a total of 36 
channels. 

Yet another embodiment of the present invention provides 
a high level of reliability in the severe radiation environment 25 
present in high Earth orbits. Specifically, the receiver 100 
(with reference to FIG. 7) withstands a total dose radiation 
level of approximately 1 00 krad. This exposure is with no box 
shielding; however, spot shielding is permissible. Moreover, 
the receiver 100 preferably is tolerant to an approximate 37 
MeV-cm 2 /mg exposure; with no single event upsets (SEUs). 
Additionally, the receiver 100 preferably is single event latch- 
up (SEL) and single event burnout (SEB) immune up to 
approximately 90 MeV-cm 2 /mg. 

Selecting parts that provide the required performance and 
survivability is a significant design challenge. The following 
are examples of the various individual components used in 
forming the embodiments of the invention. While, specific 
parts are identified, those skilled in the art would readily 4Q 
recognize that comparable parts may be used as a substitute 
for the parts identified below, and the embodiments of the 
invention are not limited to a particular type of part/compo- 
nent. One of skill in the art will also recognize that some or all 
of the hardware components may alternately be implemented 45 
in software. 

The flight radio frequency (RF) analog front end 102 (of 
FIGS. 7 and 9) may be built around a PE8510xLl/L2 GPS 
front end ASIC available from Peregrine Semiconductor Cor- 
poration, California, USA. For the FPGAs, a RTAX-2000 50 
available from Actel Corporation, California, USA, may be 
used. The combinatorial and sequential logic of Actel’ s 
device satisfactorily meets the radiation requirements of 
operation in a space environment. The random access memo- 
ries (RAMs) on the FPGA offered by Actel are actually not 55 
immune from SEUs; however, Actel offers different Error 
Detection and Correction (EDAC) algorithms in their tool set 
to increase the data resiliency of SRAM. The flight SRAM is 
preferably a 4 SRAM die available from BAE Systems North 
America, Virginia, USA, packaged in a Multi-Chip Module 60 
(MCM). 

The baseline flight oscillator is preferably an ovenized 
crystal oscillator (OXO), however a high quality temperature 
controlled crystal oscillator (TCXO) may be used for appli- 
cations that desire to trade some power, mass, and cost saving 65 
for slightly reduced performance. An RH-CF5208 ColdFire 
processor available from Motorola, Illinois, USA may be 
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used as the microprocessor 105 (of FIGS. 7 and 9). It is a true 
embedded processor with very low-power consumption and 
almost no glue logic. 

A software embodiment of the invention comprises two 
major components, as illustrated in the block diagram in FIG. 
10. The first component of the software is the low-level func- 
tions that interface directly with optional hardware compo- 
nents to produce raw code phase, carrier Doppler, and carrier 
phase measurements. These functions include the control of 
the acquisition engine 121 as well as the tracking loops 123, 
which preferably are also specialized for the weak signal 
environment. These tasks preferably operate at the interrupt 
level and have direct communication with the hardware. 

is the embodiment also preferably includes basic naviga- 
tion software, which is preferably run on the Nucleus™ real- 
time operating system available from Accelerated Technol- 
ogy, Inc. This software forms measurements, provides the 
standard position, velocity, and time point solutions when 
four or more satellites are being tracked, and handles com- 
manding and telemetry messages. Additionally, this software 
embodiment of the invention preferably includes an attitude 
determination capability when setup with a suitable antenna 
configuration. The GPS Enhanced Onboard Navigation Sys- 
tem (GEONS) 124 preferably is integrated with the receiver 
software to provide onboard orbit determination capabilities, 
and to provide accurate state estimation/propagation during 
periods of limited/no GPS observability. The navigation soft- 
ware preferably is hardware independent by utilizing a pipe- 
wall construct 125 to communicate with the low-level soft- 
ware functions (navigation 126, communications 127 (to a 
spacecraft), and the Ephermeris/Almanac 128, which pro- 
vides GPS data). As a result, the new software can be run with 
both the hardware and the receiver cards (not shown), which 
optionally allows the software and hardware embodiments of 
the invention to function in parallel. Thus, embodiments of 
the invention software optionally can be used with other types 
of hardware platforms. One of skill in the art will recognize 
that some or all of the harware functions depicted in FIG. 10 
alternately may be implemented in software. 

Yet another embodiment provides a GPS receiver 100 that 
is a fully space flight qualified GPS receiver 100 optimized 
for fast signal acquisition and weak signal tracking. The fast 
acquisition capabilities preferably provide exceptional TTFF 
with no a priori receiver state or GPS almanac information, 
even in the presence of high Doppler shifts present in LEO (or 
near perigee in highly eccentric orbits). The fast acquisition 
capability also makes it feasible to implement extended cor- 
relation intervals and therefore significantly reduce the acqui- 
sition threshold of the embodiments of the invention. 

The foregoing description of the specific embodiments will 
so fully reveal the general nature of the invention that others 
can, by applying current knowledge, readily modify and/or 
adapt for various applications such specific embodiments 
without departing from the generic concept, and, therefore, 
such adaptations and modifications should and are intended 
to be comprehended within the meaning and range of equiva- 
lents of the disclosed embodiments. It is to be understood that 
the phraseology or terminology employed herein is for the 
purpose of description and not of limitation. Therefore, while 
the embodiments of the invention have been described in 
terms of preferred embodiments, those skilled in the art will 
recognize that the embodiments of the invention can be prac- 
ticed with modification within the spirit and scope of the 
appended claims. 



21 


US 7,548,199 B2 


What is claimed is: 

1. A global positioning system (GPS) receiver comprising: 

a GPS signal acquisition component comprising a fre- 
quency domain correlation module, 

wherein said GPS signal acquisition component is adapted 5 
to acquire a GPS signal by receiving data from said GPS 
signal and processing said data to detect said GPS signal 
within one second for GPS signals with an EIRP of at 
least -160 dBW, and 

wherein said GPS receiver is adapted to operate in a space to 
environment. 

2. The GPS receiver of claim 1, further comprising: 

a plurality of GPS signal tracking correlators adapted to 
track said GPS signals with an EIRP of at least -175 
dBW. 15 

3. The GPS receiver of claim 1, wherein said GPS signal 
tracking correlators comprise a field -programmable gate 
array (FPGA). 

4. The GPS receiver of claim 1, wherein said GPS signal 
acquisition component comprises a field -programmable gate 20 
array (FPGA). 

5. The GPS receiver of claim 1, wherein said GPS signal 

acquisition component is adapted to acquire a GPS signal 
within one minute for GPS signals with an EIRP of at least 
-180 dBW. 25 

6. The GPS receiver of claim 1, wherein said GPS signal 
acquisition component is adapted to perform a Fast Fourier 
Transform (FFT) correlation process on said GPS signals. 

7. The GPS receiver of claim 1, wherein said GPS signal 
acquisition component is adapted to acquire a GPS signal 30 
without a priori data for GPS signals with an EIRP of at least 
-180 dBW. 

8. A global positioning system (GPS) receiver comprising: 

a GPS signal acquisition component comprising a fre- 
quency domain correlation module, 35 

wherein said GPS signal acquisition component is adapted 
to acquire said GPS signals with an EIRP of at or below 
-180 dBW without a priori data, and 

wherein the acquired GPS signal is acquired by receiving 
data from said GPS signal and processing said data to 40 
detect said GPS signal substantially in real time. 

9. The GPS receiver of claim 8, further comprising: 

a plurality of GPS signal tracking correlators adapted to 
track said GPS signals with an EIRP of at least -175 
dBW. 45 

10. The GPS receiver of claim 9, wherein said GPS signal 
tracking correlators comprising a field-programmable gate 
array (FPGA). 

11. The GPS receiver of claim 8, wherein said GPS signal 
acquisition component comprises a field-programmable gate 50 
array (FPGA). 

12. The GPS receiver of claim 8, wherein said GPS signal 

acquisition component is adapted to acquire a GPS signal 
within one minute for GPS signals with an EIRP of at least 
-180 dBW. 55 

13. The GPS receiver of claim 8, wherein said GPS signal 
acquisition component is adapted to acquire a GPS signal 
within one second for GPS signals with an EIRP of at least 
-160 dBW. 

14. The GPS receiver of claim 8, wherein said GPS signal 60 
acquisition component is adapted to acquire a GPS signal 
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within one minute for GPS signals with an EIRP of at least 
-180 dBW, and to acquire a GPS signal within one second for 
GPS signals with an EIRP of at least -160 dBW. 

15. The GPS receiver of claim 8, wherein said GPS signal 
acquisition component is adapted to perform a Fast Fourier 
Transform (FFT) correlation process on said GPS signals. 

16. A global positioning system (GPS) receiver compris- 
ing: 

a GPS signal acquisition component comprising a fre- 
quency domain correlation module, 

wherein said GPS signal acquisition component is adapted 
to acquire said GPS signals with an EIRP of at or below 
-160 dBW without a priori data, and 

wherein the acquired GPS signal is acquired by receiving 
data from said GPS signal and processing said data to 
detect said GPS signal substantially in real time. 

17. The GPS receiver of claim 16, further comprising: 

a plurality of GPS signal tracking correlators adapted to 
track said GPS signals with an EIRP of at least -175 
dBW. 

18. The GPS receiver of claim 17, wherein said GPS signal 
tracking correlators comprising a field-programmable gate 
array (FPGA). 

19. The GPS receiver of claim 16, wherein said GPS signal 
acquisition component comprises a field-programmable gate 
array (FPGA). 

20. The GPS receiver of claim 16, wherein said GPS signal 
acquisition component is adapted to acquire a GPS signal 
within one minute for GPS signals with an EIRP as low as 
-180 dBW. 

21. The GPS receiver of claim 16, wherein said GPS signal 
acquisition component is adapted to acquire a GPS signal 
within one second for GPS signals with an EIRP of at least 
-160 dBW. 

22. The GPS receiver of claim 16, wherein said GPS signal 
acquisition component is adapted to acquire a GPS signal 
within one minute for GPS signals with an EIRP of at least 
- 1 80 dBW, and to acquire a GPS signal within one second for 
GPS signals with an EIRP of at least -160 dBW. 

23. The GPS receiver of claim 16, wherein said GPS signal 
acquisition component is adapted to perform a Fast Fourier 
Transform (FFT) correlation process on said GPS signals. 

24 . A method of acquiring global positioning system (GPS) 
signals, said method comprising the steps: 

receiving GPS signals; 

performing a frequency domain correlation process on said 
GPS signals; 

comparing the correlated GPS signals with a predeter- 
mined threshold and acquiring a GPS signal if a corre- 
lated GPS signal exceeds the predetermined threshold; 

wherein said steps of acquiring said GPS signals are per- 
formed within one minute for GPS signals as low as 
-180 dBW. 

25. The method of claims 24, wherein said steps of acquir- 
ing said (GPS) signals are performed within one second for 
GPS signals as low as -160 dBW. 

26. The method of claim 24, fiirther comprising acquiring 
said GPS signal without a priori data. 





